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The paradoxical response, when a system produces an output — current, flux 
or rotation — in the direction opposite to an applied bias or torque, seemingly 
contradicts Newton's second law. Nevertheless, such a response, termed abso- 
lute negative mobility jll-lel], can occur in a system driven far out of thermal 
equilibrium. In earlier works, the presence of strong decoherence mechanisms 
was regarded as a prerequisite for negative mobility to appear Here we 

report an experimental and theoretical study showing absolute negative mobility 
in a fully coherent quantum regime. Using a rubidium Bose-Einstein condensate 
loaded in a time-periodically modulated optical lattice, we observe a negative re- 
sponse from a coherent quantum system. The experimental results are explained 
in terms of eigenstates of a time-periodic Hamiltonian. Our findings open new 
possibilities for the control of quantum transport in the decoherence-free limit. 

Achieving control of the dynamics of quantum systems is one of the primary goals in 
the field of nanoscience jsl- While the quest for control of nano- and microscopic systems 
is usually pursued using static confining potentials, the use of time-dependent potentials 
is less explored. Ultracold atoms provide attractive model systems for the application of 
tunable potentials [9], and experiments with atoms in time- dependent optical potentials 



have demonstrated a variety of possibilities to control atom dynamics |10l-ll3|. 

Modulated confinements also offer an entirely new perspective, allowing one to focus on 
the properties of non-equilibrium quantum systems. It is known that not only system char- 
acteristics, but also the system response to an external perturbation can drastically differ 
from those exhibited by the same system near equilibrium 1^. A counterintuitive negative 
response, namely that a system reacts in the direction opposite to the direction of applied 
force, may appear due to the interplay of three key factors. These are (i) nonequilibrium 
conditions, (ii) a nonlinearity, realized, for example, with spatially periodic potentials, and, 
as it has been reckoned to date, (iii) the presence of decoherence mechanisms . For 

instance, regimes of absolute negative mobility (ANM) have been detected in solid-state 



electronic devices 



^ and have also been investigated theoretically with different classical 
stochastic models [4,5]. Recently, absolute negative mobility effects have been observed in 
situ, by tracking an ensemble of charged polystyrene beads with a real-time video micro- 
scope (6|. 

Here we show that a fully coherent single-particle quantum dynamics can exhibit negative 
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response, which persists on experimentally relevant timescales. We present a theoretical 
analysis via the Floquet formalism [l5| and outline the mechanisms which give rise to 
ANM quantum regimes. The theoretical predictions are corroborated with the experimental 
results obtained by measuring the velocity of dilute rubidium atomic BEC loaded into a 
time-modulated periodic optical potential and subjected to a net bias. 

To start, consider the model of a single particle of mass M placed into a time- and 
spatially periodic potential, V{x)A{t) with V{x + L) = V{x), A(t + T) = A(t), of periods L 
and T respectively. In addition, assume that the particle is exposed to a tunable net bias, 
F. The corresponding Hamiltonian reads 

H{x,p,t)=p'^/2M + V{x)A{t) + Fx. (1) 

The undriven system obeys A{t) = 1 and F = 0; the Hamiltonian then is spatially- 
periodic, and the system quasimomentum space is spanned by the Bloch bands, 
K G [—7t/L,7i/L] and n = 1,2,.... By assuming that the initial state is well-localized at 
the point k = 0, and considering a weak bias F as a force which pulls the system through 
the quasimomentum space, we arrive at a linear time- dependence of the quasimomentum 
K{t) = Ft/h. The response of the system to a weak bias on a time scale t < 2TTh/FL = Tb, 
where Tg is the period of Bloch oscillations, is determined by the curvatures of the bands, 
h~'^d'^En{k) / dn^ , i. e. by the effective mass. Hence a negative mobility, following the recipe 



given by Kromer 16|, demands a negative effective mass. This implies that an ANM re- 
sponse from the ground state of a one-dimensional lattice is rendered impossible, because 
the curvature is always positive near the center of the first Brillouin zone, see Fig. la. 

The dynamics of time-modulated systems is more diverse. The system Hamiltonian is 
periodic in time, and the solution of the corresponding Schrodinger equation for a given k 
can be obtained as a set of eigenf unctions, lipn^nii))- The eigenf unctions satisfy the Floquet 
theorem, so that iV^n^W) = exp(-2e„(K)t//i) |0„,«;(O)), where \(j)n,K{t + T)) = \4>n,^,(t + T)) . 
The quasienergies e„(/t) are conventionally restricted to the interval [—tnr/Tjtnr/T]. The 
right upper part of Fig.[T]D depicts a typical Floquet band structure of a driven optical lattice. 
It exhibits a complex, web-like structure, which is much more tangled than the Bloch band 
structure of the undriven system. The set of Floquet states can be ordered in different ways. 
In the present context it is advantageous to order them with respect to their mean energy 
at the point k = 0, En = {{(pnfl\H\(pnfl))T [LL\- The Floquet bands are well separated. 
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FIG. 1: Responses to a weak bias of a quantum particle placed into an undriven and 
ac-driven periodic potentials, (a) Typical Bloch band structure of an undriven system, with 
the ground state marked by a dashed line. The action of a weak bias can be represented as a 
slow motion through the quasimomentum space. The velocity of a particle localized at a certain 
point in the quasimomentum space is proportional to the slope of the band at the localization 
point. A particle, initially located in the ground state at the point k = 0, is capable of producing 
a weak positive response only. Like a ball placed on a hillside, it always rolls downhill, (b) The 
Floquet spectrum of the driven system. The average velocity is proportional to the local slope of 
the corresponding quasienergy band. The thick red line marks a ballistic diabatic band, which is 
the state of lowest mean energy at the point k = 0. A quantum particle, initially prepared in this 
state, produces a strong negative response when subjected to a weak bias, i.e. its center of mass 
moves into the direction opposite to the applied bias. This effect can be thought as the uphill 
motion of a bouncing ball. 

and, although they may approach each other at the points of avoided crossing Q, their 
order is preserved over the entire Brillouin zone. An adiabatically slow motion through 
the quasimomentum space typically respects all avoided crossings. However, for any finite 
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force, F, there exists always a set of narrow avoided crossings such that the system ignores 
them when moving across. The corresponding trajectories in /t-space form so-called diabatic 
bands [itI. Once populated, a diabatic state would remain decoupled from other Floquet 
states during the excursion through the quasimomentum space. 

The transport properties of Floquet states are characterized by their average velocities, 
Vn{i^) = {vn.K.{ t))T , which cau be calculated from the corresponding quasienergy, Vn{f^) = 
hr^den{K) / dn 18|. When the system Hamiltonian is invariant under time reversal, t — t- — t, 
all Floquet bands are flat at the center of the Brillouin zone, thus assuming zero average 



velocities 



19| . Outside the point k = the situation is different: Even in the case when 



the time-reversal symmetry is present, the Floquet states acquire nonzero velocities upon 
departure from the center, and some of them form ballistic diabatic bands. 

Consider next the situation when at time t = one starts an experiment with a dilute 
cloud of ultracold atoms, smeared over many wells of an optical potential. An initial wave 
function in form of a wave-packet localized at the point /t = will overlap predominantly with 
the Floquet ground state of a driven lattice, i. e., with the Floquet state with lowest energy. 
El. Different from the Bloch ground state, the Floquet ground state may itself represent a 
ballistic diabatic band outside the center of the Brillouin zone, a situation shown in Fig. lb. 
Alternatively, a Floquet ground state can be set in contact with one of the ballistic bands, 
through an avoided crossing at the vicinity of /t = 0. Then, the corresponding ballistic band 
will also be populated, and will determine the system mobility response, provided that the 
band was populated substantially and (or) the band velocity is large enough. If, in addition, 
the band velocity is negative, one will obtain an ANM response. 

A general recipe to observe an ANM response is the following: Firstly, the system must 
be in a deep quantum regime, ideally with only a few Floquet states within the potential 
range. Otherwise, the atomic cloud would be distributed among many eigenstates, and the 
system response will consist of many uncorrelated contributions. Secondly, the modulation 
frequency u should be close to the resonance frequency, ur = Air'^h / L"^ M . The resonance 
modulation will drive the system deeply into a nonequilibrium regime, where the structure 



of the Floquet spectrum is maximally tangled |17l . Il9l | . The frequency region around the 
resonance frequency uji should then be explored systematically in order to locate ANM 
regions. There is a tip-off that helps in this search. Namely, once one encounters a strong 
positive response, one can be confident that there is a strong negative response located 
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FIG. 2: Negative mobility of rubidium atoms, (a) Time of flight images showing the atomic 
momentum distributions after t = 26T periods of modulation for F = —0.0181 E^/X (left) and 
F = 0.0181 Er/\ (right), which is near the maximal observed ANM response respectively, and in 
the absence of a bias force (center). The used modulation frequency is ijo/ojr = 1.030. The visible 
atomic diffraction orders are s = —2,..., 2 (from bottom to top), (b) Measured mean velocity of 
the atomic cloud, in units of hk/Mjib, versus applied force F, in units of Er/X. The error bars 
show the standard deviation of the mean. The blue data points correspond to the data shown in 
(a). The solid lines depict the results of a numerical simulation of the experiment for an initial 
state in the form of a narrow Gaussian packet in the quasimomentum space, with the center at 
K = and width of fj^ = 0.04^/c. 

nearby. This is so because the ballistic states with opposite velocities are always paired in 
systems obeying time-reversal invariance. 

In our experiments we use a time-modulated optical potential, U{x,t) = V{x)A(t), 

V{x) = Vi/2cos{2kx) + V2/2sm{4:kx), (2) 
A{t) = Ai sin2(wt/2) + A2 cos^(a;t) , (3) 

where A is the wavelength of the driving laser field and k = 2tt/X. With the spatial period- 
icity L = A/2, we arrive at a resonance frequency of ur = Sur, where Ur = hk'^/2Mfib is the 
recoil frequency. For the used ^''Rb atoms and A = 783nm, we have Ur = 27i ■ 3.71kHz. The 
parameters used in experiments are: Vi = S.IE^, V2 = 0.7Er, where = TiUr is the recoil 



energy, and Ai = 0.7, A2 = 0.44. Compared to a standing wave lattice, our biharmonic 
lattice provides an enhanced control of the system dynamics, while leaving the system per- 
fectly time-symmetric. A rubidium Bose-Einstein condensate is initially freely expanded for 
2.5ms, after which its interaction energy has transformed into kinetic energy. Subsequent 
velocity selection reduces the momentum spread to a value of ±0.1hk, and the ballistically 
falling atoms are exposed to a temporally varying optical lattice potential with an additional 
bias force. During the free fall of the atoms, one of the lattice beams is acousto-optically de- 
tuned with a constant chirp-rate to accelerate the lattice with respect to the lab frame 20 |. 
In this way, we can accelerate the lattice to leave it stationary with the ballistically down- 
wards falling atoms, or leave a variable residual accreleration that acts like an external force. 
After the interaction with the driven lattice potential of variable net bias, the atoms freely 
expand for 15ms, and an absorption image is recorded. With this time-of-flight technique, 
the atomic velocity distribution is analyzed. 

Typical experimental data, as presented in Fig. 2a, consist of several atomic diffrac- 
tion peaks. The mean velocity of the atomic cloud is calculated as v = p/Mnb with 
p = 2hkY^^s\cs\'^, where |csp denotes the fraction of atoms in the s-th order momentum 
state, \2shk), with s = ±1,±2. The parameters used here were a modulation frequency 
oj = l.OSuji, an interaction time corresponding to t = 26T periods of the modulation and a 
net acceleration of F = —0.4,0, and OAMptb ■ m/s'^, which correspond to the dimensionless 
values FX/ Er = —0.0181, and 0.0181 for the left, middle, and right image of Fig. 2a respec- 
tively. The measured mean velocity of the atomic cloud versus the external force F in units 
of Er/X is shown in Fig. 2b. For small values of F, negative atomic mobility, corresponding 
to a motion against the external force, is clearly visible. The two data points marked in 
blue correspond to the data with F = ±0.0181 Er/X, see the diffraction data of Fig. 2a, 
left and right, respectively. In other measurements, we have studied the atomic transport 
versus the drive frequency u at fixed external force. Corresponding data is shown in Fig. 3a 
for F = —0.0181 Er/X, where absolute negative atomic mobility is observed for values of 
uj/ur roughly between 1.02 and 1.15. For smaller frequencies a strongly positive mobility is 
observed, in accordance with the theoretical prediction that ballistic bands always occur in 
pairs. 

An important issue is the stability of the ANM response. In this respect, it is noteworthy 
to confront the phenomenon at two different limits, classical and quantum ones. In the 
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FIG. 3: Control of mobility, (a) Mean atomic velocity as a function of the modulation frequency 
(J for a fixed bias force, F = —0.0181 Er/X, and exposition time 26T. (b) Time evolution of the 
mean velocity for two frequencies of the driving, oj/ujr = 1.05 (red dots) and oj/ujr = 1.076 (blue 
squares). The solid lines depict the results of a numerical simulation of the experiment for an 
initial state in the form of a narrow Gaussian packet in the quasimomentum space, with the center 
at K = and width of = O-OAhk. 



classical dissipationless limit, a stationary motion against constant bias is possible j2l| due 
to the existence of invariant manifolds, i.e. transporting regular islands, periodic orbits, and 
cantori, in the phase space of an ac-driven Hamiltonian system 22|. Diabatic bands can be 
viewed as the quantum counterpart of classical ballistic manifolds. However, this analogy 
is not exact. In contrast to classical manifolds, diabatic bands are not completely isolated, 
since any finite bias sets the band into contact with other states. All avoided crossings have 
finite widths, and the Landau-Zener tunneling between Floquet states (l5| will eventually 
lead to the regime of acceleration along the applied bias. Nevertheless, the relaxation to 
this asymptotic regime can be extremely slow, thus transforming the negative response into 
a long-lasting metastable phenomenon. In order to confirm this prediction, we measured 
the velocity of the atomic cloud as a function of the exposition time for two different values 
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of the driving frequency. Fig. 3b presents corresponding results of the measurements. For 
a modulation frequency oj/ujr = 1.05 (red dots), the mean momentum of the atomic cloud 
increases initially, reaching a broad maximum value of OAhk after around 30 modulation 
periods, and then starts to decrease again, approaching zero at t ~ SOT. When slightly 
changing the modulation frequency to u/ur = 1.076 (blue squares), the observed atomic 
momentum increases to roughly 0.9hk, with this value of the atomic transport beeing main- 
tained to much longer modulation times. This finding is in agreement with the fact that 
the structure of Floquet bands is very sensitive to the system parameters, and the atomic 
response can be controlled over a wide range. 



Our results can be interpreted as a realization of a quantum engine [23|, in which the 
system self-adjusts itself to perform a work against an external force. The observed quantum 
ANM effect can not only be seen as a next step towards fully controllable transport of 
ultracold matter {2^, 25I, but also highlights the salient fact that single-particle quantum 
systems can produce a spectrum of nontrivial phenomena, when moved far from equilibrium. 
For the future, we expect that inclusion of many-body effects of interacting condensate 
atoms [9I, will enrich the present spectrum even further. 

This work was supported by the DFG Grants Wel748/7 (M.W.), HA1517/31-2 (S.D. and 
P.H), the German Excellence Initiative "Nanosystems Initiative Munich (NIM)" (A.V.P., 
S.D. and P.H.). 

Supplementary Material 

Experimental set-up and procedure. The experimental set-up was described previ- 



ously in Refs. 
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20| . A spin-polarized Bose- Einstein condensate with about 6 ■ 10^ ^^Rb 
atoms in the Zeeman sub level m^r = is created by evaporative cooling in a quasistatic dipole 
trap with an applied small magnetic quadrupole field. The dipole trap is then switched off 
and the cloud performs a ballistic free fall in the earth's gravity field for 2.5ms, where in- 
teraction energy is transformed into kinetic energy. The condensate is then exposed to a 
350 /is-long Bragg pulse, which transfers a narrow slice of atoms with momentum spread 
±0.1hk into a momentum state \2hk), followed by a 150 /is-long Raman pulse transferring 
the momentum- reduced sample into the zero momentum state \Ohk) with mp = —1. The 
velocity selected atoms are subsequently subjected into a vertical Fourier-synthesized optical 
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lattice potential that is periodically modulated in time. The lattice potential is generated 
by superimposing two lattice harmonics with spatial periodicities A/2 and A/4, with the 
latter being generated by a multiphoton Raman technique 20[ . Here A = 783nm is the laser 
wavelength, some 3nm red-detuned from the rubidium D2-line. The experiments are done 
in ballistic free fall with the gravity force in the atomic frame partially compensated in a 
controllable manner by linearly chirping the frequency of one of the laser beams. After a 
fixed exposition time, all laser fields are switched off and the atoms freely expand for 15ms, 
after which an absorption image is recorded onto a CCD camera. 

Numerical methods. We model the dynamics of a wave-packet, initially localized at 
the center of the BriouUin zone, with a wave function |\E'[k1) = -7= — e~'*^/^°"«(iK. This i 

' I L J/ ^/27^<T.. J— 00 



IS 



a good approximation of a momentum slice used in the experiments, and the dispersion 
serves as a fitting parameter. The initial wave-packet has been sampled in a first Brillouin 
zone with 5000 points. Every slice is then propagated by using the integration scheme from 



Ref. 



19|. Finally, the contribution from all points are summed up. Numerical parameters 



correspond to the values given in the main text. 
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